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LONG-TERM  GOALS 

A  better  understanding  of  oceanic  variability  via  modeling  studies  of  circulation,  entrainment,  mixing 
and  convection  in  the  coastal  ocean.  Development  and  use  of  high-resolution  models  for  the  study  of 
dynamic  processes  and  for  the  investigation  of  specific  oceanic  regions. 

OBJECTIVES 

The  particular  objective  of  this  project  is  to  understand  the  physics  of  the  mesoscale  motions  across  the 
Adriatic  Sea  well  enough  to  simulate  them  faithfully  in  high-resolution  models. 

APPROACH 

The  selected  model  is  DieCAST  because  of  its  extremely  low  level  of  dissipation  at  the  grid  scale 
(about  2  km  in  the  horizontal),  and  forcing  fields  are  interpolated  from  COAMPS  data  sets  (from  the 
NRL  in  Monterey).  In-situ  data  used  for  validation  are  surface  drifter  trajectories  (from  OGS  in 
Trieste).  Model  results  are  sampled  where  drifter  data  exist,  and  identical  statistics  are  performed  on 
both  drifter  speed  and  simulated  currents.  In  particular,  mean  and  turbulent  kinetic  energies  are 
calculated  and  compared.  The  project  is  integrated  in  the  ONR-sponsored  Dynamics  Of  Localized 
Currents  and  Eddy  Variability  In  The  Adriatic  (DOLCE  VITA)  Program. 

WORK  COMPLETED 

The  DieCAST  ocean  model,  applied  in  a  high-resolution  mode  to  the  entire  Adriatic  Sea,  was  used  to 
simulate  surface  currents  during  periods  for  which  there  is  better- than-average  drifter  trajectory 
coverage,  so  that  drifter  and  model  statistics  could  be  compared.  These  periods  were:  1-23  October 
2002,  1-25  February  2003,  and  1-30  June  2003.  Interestingly  enough,  some  of  the  selected  periods 
coincided  with  strong  wind  events,  either  bora  and  sirocco.  The  corresponding  COAMPS  datasets 
were  obtained  and  mapped  onto  the  DieCAST  model  grid.  Prior  model  spin-up  was  performed  with 
climatological  fields  for  3  months.  Forcing  was  then  switched  to  the  COAMPS  data,  and  the 
simulation  continued  for  a  month  or  more.  During  the  selected  comparison  periods,  surface  currents 
were  sampled  from  the  model  results  at  times  and  places  were  drifter  velocities  are  known  so  as  to 
create  two  surface  velocity  sets  of  identical  size.  At  each  available  location  (5x5  model  grid  box  with 
5  or  more  drifter  data),  means  and  variances  were  calculated  from  the  time  series  and  transformed  into 
mean  kinetic  energy  (MKE)  and  turbulent  kinetic  energy  (TKE),  per  mass  (units  of  cm“/s  ). 
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RESULTS 


Comparisons  of  MKE  and  TKE  between  model  simulations  and  drifter  trajectories  are  favorable,  but, 
of  course,  one  looks  for  differences  more  than  similarities  in  order  to  obtain  insight  about  ways  to 
improve  the  model.  For  October  2002,  a  period  during  which  the  Western  Adriatic  Current  (WAC) 
was  particularly  strong,  drifter  data  show  high  values  of  both  MKE  and  TKE  all  along  the  Italian 
Coast,  and  the  model  behaves  likewise  except  that  its  MKE  is  lower  and  TKE  higher  (Figure  1).  In 
order  words,  the  model  overestimates  the  WAC’s  meandering  and  eddying.  This  indicates  that  the 
model  is  insufficiently  dissipative. 

For  February  2003,  a  time  of  strong  bora  wind,  the  model  performs  nearly  identically  to  the  drifter 
data,  except  at  the  southern  tip  of  Istria  and  along  a  line  extending  southwestward  from  there.  This  is 
where  the  bora  generated  an  intense  cross-cutting  jet,  and  the  model  underestimated  its  strength.  The 
model  is  not  wrong  per  se,  but  its  “surface”  velocity  is  by  design  the  velocity  averaged  over  the  top  5m 
of  the  water  whereas  drifters  sample  the  velocity  in  the  only  top  lm  of  the  water.  Naturally,  a  swift 
wind  creates  a  drift  current  highly  sheared  in  the  vertical,  and  the  difference  between  lm  and  5m  can 
be  significant.  The  remedy  is  to  increase  the  model’s  vertical  resolution  near  the  surface  so  that  its  top 
velocity  is  more  representative  of  what  drifters  actually  measure. 

For  June  2003,  a  time  of  pronounced  cross-basin  jets  and  gyres,  the  comparison  between  model  results 
and  drifter  data  is  excellent,  with  a  slight  degree  of  overestimation  of  TKE  by  the  model,  suggesting 
once  again  that  the  model  is  insufficiently  dissipative. 

IMPACT/APPLICATIONS 

The  results  demonstrate  that  effective  simulations  of  the  mesoscale  variability  of  the  Adriatic  Sea 
require  a  low-dissipation  model,  a  grid  resolution  of  at  least  2  km,  realistic  surface  forcings  (wind 
stress,  heat  flux,  etc.)  of  similar  spatial  resolution,  AND  a  turbulence  closure  for  near-surface  mixing. 
This  ought  to  impact  the  Mediterranean  Sea  models  currently  used  by  the  US  Navy.  It  also  leads  us  to 
anticipate  success  should  DieCAST  be  later  configured  with  data  assimilation. 

RELATED  PROJECTS 

This  project  was  a  component  of  the  multi-project  DOLCE  VITA  Program  funded  by  ONR-PO  and 
focusing  on  the  Adriatic  Sea.  Related  projects  are  those  of  Craig  Lee  (Univ.  Washington,  TriSoarus 
towed  profiling),  Pierre-Marie  Poulain  and  Elena  Mauri  (OGS-Trieste,  surface  drifters),  and  Mirko 
Orlic  (Univ.  Zagreb,  East  Adriatic  Coastal  Experiment). 
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Figure  1.  Comparison  of  surface  mean  kinetic  energy  (MKE  -  top  row  -  values  ranging  from  about 
zero  to  300  cm2/s2)  and  turbulent  kinetic  energy  (TKE  -  bottom  row  -  values  ranging  from  about  10 
to  100  cm2 Vs2)  calculated  from  both  drifter  trajectories  (left  panels)  and  DieCAST  model  simulations 
(right  panels)  for  the  period  1-23  October  2002.  Blank  areas  indicate  regions  of  insufficient  drifter 
data  (less  than  5  per  box).  This  period  was  one  of  intense  flow  down  the  Italian  coastline  and  of 
complex  flow  in  the  northern  basin.  The  model  faithfully  reproduces  regions  of  strong/weak  mean 
flow  and  low/high  variability.  However,  discrepancies  in  values  (MKE  lower  in  model  by  10%  and 
TKE  higher  in  model  by  20%)  indicate  that  the  model  overestimates  the  meandering  and  eddying  of 
the  coastal  current  along  Italy,  and  this  is  attributed  to  insufficient  dissipation 

on  the  part  of  the  model. 
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